't Hooft and Wilson loop ratios in the QCD plasma 



1 Introduction 

The plasma phase of QCD is nowadays object of experimental study at RHIC. 
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i The spatial 't Hooft loop measuring the electric flux and the spatial Wilsonloop 

■ measuring the magnetic flux are analyzed in hot SU(N) gauge theory. Both display 

area laws. On one hand the tension of the 't Hooft loop is perturbatively calculable, 
» I ' in the same sense as the pressure. We show that the 0(g 3 ) contribution is absent. 

f"S |. The ratio of multi-charged 't Hooft loops have a remarkably simple dependence 

, on the charge, true up to, but not including, 0(g i ). This dependence follows 

> ■ also from a simple model of free screened colour charges. On the other hand the 

' surface tension of the Wilsonloop is non-perturbative. But in a model of screened 

free monopoles at very high temperature the known area law follows. The density 
of these monopoles starts to contribute to 0(g 6 ) to the pressure. The ratio of the 
i multicharged Wilson loops is calculable and identical to that of the 't Hooft loops. 

> 
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^ The equilibrium properties of the plasma have been analyzed by numerical simulations and 

analytic studies. The latter mainly through the use of perturbation theory. The first obvious 
question is how we can distinguish between a plasma phase and a hadron phase. 

One order parameter is the free energy of a single quark or Wilson line. Another the free 
energy of a single Dirac monopole. The first one is the order parameter for a discrete global 
symmetry, called electric Z(N) symmetry. This symmetry is a symmetry of the Euclidean 
path integral, as long as only Z(N) neutral fields appear. It is broken in the high T phase, as 
signalled by the non vanishing VEV of the Wilson line. The symmetry of the Dirac monopole 
line is always broken. 

There are also two canonical order parameters in QCD : the spatial Wilson loop, and the 
spatial 't Hooft loop. The first obeys an area law as numerical studies have shown @, both in 
the confined and in the plasma phase. This corresponds-ta the Dirac monopole loop having a 
VEV in both phases. This is in accord with the old ideaBE9 that there is a screening monopole 
condensate. 

The 't Hooft loop has perimeter behaviour in the cold phase (if quarks are absent), and 
area behaviour in the plasma phase. So the 't Hooft loop is the canonical order parameter 
that distinguishes the two phases. This corresponds to the free energy of a heavy quark being 
infinite in thje cold and finite in_the plasma phase. In the latter phase it has been calculated 
analyticallyH and by simulationcl. 

What sets the 't Hooft loop and the Wilsonloop apart physically, is that the 't Hooft loop 
measures the electric colour flux in the plasma, whilst the Wilson loop measures the magnetic 
colour flux. The plasma is par excellence containing free or almost free screened colour charges. 
In such surroundings the_aijea law exp— a(T)A of the 't Hooft loop follows by a well-known 
simple physical argument E3'B to be: 

a(T) ~ l E n E (T) (1) 

with the Debye screening length I e multiplying the density of electric colour charge due to the 
gluons. All what goes into the area law is that gluons are screened, free and obey a simple 
thermodynamical distribution law. 
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Though at extremely high temperatures (extreme is here 10 3 GeV, so rather academic) 
simple minded perturbation theory works for the free energy, it does not work at temperatures 
below. Resummation techniques have been devised to repair the situationcl. 

However to work for the Debye mass one has to go to even higher temperatures. So the 
situation is different from one observable to another. Though the reason for the bad convergence 
is not well understood, the facts are the following: the perturbative contributions coming 
exclusively from the excitations of order T do converge well. The trouble with the convergence 
starts when taking into account contributions from scales comparable to the electric screening 
mass O(gT) and the magnetic screening mass 0(g 2 T). For the free energy -for example the 
0(g 3 ) renders the convergence quite bad. For a systematic view on this, see ref.EZI and references 
therein. n 

One example where perturbation theory may work better is that of the spatial 't Hooftlil 
loop, because there 0(g 3 ) are absent as we will show in this paper. So we may hope that our 
ratios are not far off its value at the phase transition, and can play a role in the determination 
of the universality class of the 3d spin model at T c . 

In this paper we discuss in the first place some quantitative properties of the electric surface 
tension in the two loop approximation. In particular we show that the surface tension of the 't 
Hooft loop is identical to the surface tensioruDf Z(N) domain walls. The latter are formulated 
in terms of the potential for the Wilson lineEl. 

As is well known, the minima of the potential for the Wilson line are precisely where the 
Wilsonline takes on centergroup values. If we denote by a k the surface tension one gets by 
tunneling from the value 1 to the value expik^- through the straight path connecting them, 
then we find to one loop order and two loop order: 

Uk_ = N-k 

kai N - 1 1 ' 

We will show that there is a one to one correspondence between the a% and the surface 
tension of an 't Hooft loop with multiple center group charge k. 

So eq.(||) tells that the 't Hooft loop with the strength of k charges, has less energy than k 
"fundamental" , or singly charged 't Hooft loops. Thus fundamental 't Hooft loops will attract 
each other. The effect disappears for large N and finite k as expected in a theory without 
interactions, i.e. without correlations for loops. For k on the order of N this is not true any 
longer. 

By a trivial rewriting of eq.(^J): 

a k _ k(N-k) 

V x ~ N-l [S) 

we see that charge conjugation is respected, as obviously should be the case. This means that 
only from N=4 on the ratio can take non-trivial values. 

When we define the coupling such that the one loop effects disappear in the renormalization 
of the gauge coupling in the reduced model we have: 

(T \ = ^ff-^l - (15.2785.. - ll/3( 7 , + 1/22))^) (4) 

There is a slight modification in comparison with ref.i due to a different definition of the gauge 
coupling there. 7, = 0.577215... is Euler's constant. 
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Right at the critical temperature there is the idea of the universality class of Z(N) spin 
systems being the same as that of the SU(N) gauge theory in one dimension higher. 

For N=2 or 3 the Z(N) spin system in question is unambiguously defined, and lattice data 
have given strong support to this universality. In particular the SU(2) surface tension has the 
same critical exponent as its analogue, the order-order interface, in 3d Z(2) Ising spin systemsu. 

But for N larger than 3 the space of couplings defines more than one model and, what 
is worse, the critical behaviour of these models is not unique. For example the family of 
Potts models has weights such that for equal spins the energy is 1, whereas for unequal spins 
the energy is 0. One draws immediately the conclusion that for all order-order interfaces the 
tension is the same. So our ratios, if still valid at the transition, would exclude Potts models. 

Their transition is known to be first order. For clock models it is second order in d=3, and 
infinite order in d=2, for N > 4. For Villain models the transition is infinite order for TV > 4 
in d=2, but for d=3 we have not been able to find any litterature. 

On the other hand a recent paper E3 sees a first order transition for the SU(4) pure gauge 
theory. 

The lay-out of this paper is as follows. In section ^ we give the simple intuitive argument 
that leads to the area law for the 't Hooft loop, eq.([j]). 

Section [| gives the quantitative formalism, and explains the identity of Z(N) domain wall 
tension and surface tension of the 't Hooft loop. 

Section || analyzes the symmetries of the potential term in the effective potential and derives 
the one loop result for the ratio. 

Section pi derives the two loop result for the ratio. 

Section |6| treats the kinetic term in the effective potential in more detail to justify the 
procedure used in the section before. It is a little technical, and can be skipped by the reader 
as the rest of the paper can be read independently. We show that at two loop order the wings 
of the profile of the " soliton" spanned by the 't Hooft loop include the same non-perturbative 
effects as the Debye mass. Nethertheless the integrated profile, i.e. the surface tension, is 
finite to two loop order. It ignores to this order contributions from the electric and magnetic 
screening scales. This is like in the free energy. A new result is that we show why there are no 
0(g 3 ) corrections to the tension, unlike for the free energy. 

In section [7] we show that the straight path between minima is certainly a local minimum. 
For SU(4) numerical analysis shows it is also a global minimum. 

In section || we revisit the simple method of section || and extend it to incorporate the 
counting rules giving rise to the ratios and rederive the latter. 

In section [)] we analyze the Wilsonloop in terms of a speculative monopole model, which 
has the virtue to predict the same ratios as for the 't Hooft loops. 

In the last section we draw conclusions and view the prospects. 

To avoid clutter in the formulas we have put the coupling g = 1 . In sections |[ [| and || 
we have put the coupling back in the equations. 

2 Elementary facts 

Spatial Wilson loops W(L) have been measured since long in simulations. They are expessed 
in terms of the vector potential as: 

W(L) = TrVcxpi J dl.A (5) 
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Here A = A a \ a , the Gell-Mann generators being normalized to |. The loop, being in the 
fundamental representation, carries the Z(N) charge k = 1. 

We will be interested in Wilsonloops with multiple Z(N) charge k. That is, if k quarks are 
involved in the representation it has charge k and such a loop will-be written as Wk(L). 

Spatial 't Hooft loops Vi(L) have only recently been advocated! as useful orderparameters. 
The 't Hooft loop is defined as a gauge transformation that has a jump on a surface with 
L as border, and the jump expilj^ is in the centergroup of the gauge group SU(N). More 
cxplicitely, if Qi (x) = exp iui (x) is such a gauge transformation in the defining representation 
of SU(N), then 

Vi(L) = expiDui.E (6) 
is such an operator. Our notation is E = E a X a . The canonical commutation relations read: 

[E°(x),A b k (y)] = ±6 a > b 6 j , k 5(x-g) (7) 

This definition suffices for the 't Hooft operator to have a uniquely determined action on a 
physical state as shown in the next subsection. The spatial Wilson loop is of course explicitely 
gauge invariant. 

They obey the 't Hooft commutation relations El: 

2-7T 

V;(L)H / A;(L')V;(L) _1 = exp (ikl—n(L, L'))Wk(L') (8) 

The looping number of the two loops is written as n(L, L'). It can take all integer values. 
As is well known this commutation relation imposes that at least one of the loops must follow 
an area law El. 

The 't Hooft loop captures electric colour flux, and the Wilson loop magnetic colour flux. 
We will see in this section and in section [)] more precisely how. 

In the subsection we show how the simple formula eq. ([!]) comes about. 



2.1 A simple analysis of the 't Hooft loop tension 

The definition of the spacelike singly charged 't Hooft loop V\(L) is that of a discontinuous 
gauge transformation. We can rewrite it as a loop measuring the flux of colour hypercharge: 

47r r - - 

Vi(L) = expi— TrYE.dS (9) 

N Js(L) 

The canonical electric fieldoperator E is projected onto the hypercharge matrix 
r = diag(l,-..,l,-JV+l). 

This operator, despite the fact that it transforms under gauge transformations as an^djoint, 
has a unique and gauge invariant action on physical states, as was expounded in ref. c3. It is 
worthwhile to recall this. 

First note that the 't Hooft commutation relation is reproduced by our "flux" definition in 
eq. ([)]) , with 1 = 1. This follows from the canonical commutation relations between electric field 
strength and vector potentials. Second, physical states can contain non-local gauge invariant 
quantities in the form of spatial Wilson loops. The action of V\(L) on such a loop is fixed by 
the 't Hooft commutation relation and the looping factor. So the action on a physical state is 
the same for both expressions for the loop. 
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From this it follows there is nothing sacred about the hypercharge direction Y. A regular 
gauge transform of V\ (L) : 



U{n)V x {L)U\n) =ey^>i^ J TrStfYQE.dS (10) 

has the same effect in the 't Hooft commutation relation, so acts the same on a physical state. 
What is sacred is the normalization in front of the surface integral, because that determines 
the centergoup transformation expi^-HY^ , which equals expi 2 ^ whatever 57 is. 

The thermal average < Vi(L) > has an area law behaviour in the deconfined phase and 
it is qualitatively and quantitatively understood why □. It is simply the consequence of the 
gluonic colour charges being free, and screened. To a first approximation screening means that 
only charges can contribute to the flux, that live in a slab containing the loop of thickness Ie, 
the screening length |] The flux from one single hypercharge TV in the slab will be it because 
only one half of the flux goes through the loop L. So the value of the loop from the single 
charge equals expz7r = — 1. To understand this in a slightly different way consider the adjoint 
Wilsonloop, given by TrPexpi J L , dl.A ac ij. 

When it loops the 't Hooft loop once nothing happens because it does not feel the center- 
group, and the commutation relation becomes, using eq.(^) and the canonical commutator of 
E and A: 

V k {L)W ad3 {L')V k {L)- 1 = exp (ik^-Y adj )W adj (L') (11) 

Here the matrix Y ad j is the adjoint representation of the hypercharge Y defined underneath 
eq. (|^) . So its iV 2 — 1 diagonal matrix elements consist of the differences of the diagonal elements 
of Y, so are or ±iV, so the phase factor is always one. 

The adjoint Wilson loop sends a flux iV through the 't Hooft loop. That is why in the 
commutation relation it will just produce a factor expi^-N = 1. But the point like gluon 
charge sends only a flux through the 't Hooft loop, hence the factor —1. We will use this 
argument in more detail in section || 

Since the gluons are free the probability P{1) of 1 gluons being in the slab will determine 
the average of the loop Q 

<V 1 (L)>=^2(-l) l P(l) (12) 
i 

The probability P(l) follows from the grand canonical Gibbs distribution. Far away from a 
critical point the distribution is centered around the average I with a finite width w proportional 
to VI. For the ideal gas and the high temperature Bose gas it is consistent with a Poisson 
distribution where w = VI. Hence we will take the Poisson distribution for simplicity, keeping 
in mind the caveat of the width. 

So with the probability taken Poissonian 

P(0 = i(0'exp-r (13) 

one gets for the average of the loop 

< Vi(L) >=exp-2Z" (14) 



"This is an overestimate. See ref. for better. 

''This is a well known quantity called the probability generating function 11(2) = z l P(l) with z = — 1. 
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by combining cq.(|12[) and eq. ([L3|). I is the average number of charges in the slab. Now that 
number is the volume Ie times the area A of the loop multiplying the density n of gluons and 
antigluons that have non-zero hypercolour charge ±N: 

r~ (N- l)T 3 l E A 

Together with the known screening length Ie = — we find for the tension u\ of a 

singly charged t Hooft loop a\ ~ \y J N • 

More generally there is a temperature dependent proportionality constant v{T) in the 
exponent of eq. (jl4|) due to the distribution not being Poissonian. This constant will not affect 
the area behaviour, and it will drop out in ratios of multiply charged loops. 

Suppose we want to know what the tension oi of an 't Hooft loop with charge 2 is. Let's 
consider to that end two unit charge 't Hooft loops of equal size and parallel at a distance d. 
We calculate their correlation by the same simple method as above. If d » Ie there will be 
no correlation. However for d ~ Ie a charge in between the two loops will give a — sign to 
both loops so the correlation between the loops will have a plus sign from the whole volume in 
between. So the effective volume that does contribute will be only on the outside of the two 
loops and hence smaller than the total effective volume if the two loops are far apart. So the 
surface tension of the correlation will become smaller, the loops start to attract. 

One of the aims of this paper is to do an analytic calculation of these effects for the 't Hooft 
loop, and the outcome is qualitatively the same as. in eq.(Q). Quantitatively the corresponding 
surface tension o\ has been calculated long ago u in perturbation theory, to two loop order 

0(g 2 )- 

3 The quantitative approach 

The result for the surface tension of the 't Hooft loop is reminiscent of semiclassical physics. In 
fact it can be computed semiclassically in a systematic way. The formalism will be explained 
in this section. We will, in doing so, also show the relation with the so-called Z(N) domain 
walls B. 

It will turn out that both the profile of a domain wall, and the profile that describes the 
surface tension are given by the path ordered Wilsonline 

rP n 
P(A (x)) = Pexpz / dTA Q (x,r)= lim TT exp iAT k A (x, r k ) (15) 
J0 fc=l 

which we can diagonalize to get the gauge invariant phases of the loop. As is well known its 
Euclidean average can be rewritten as the free energy of a heavy quark in the fundamental 
representation |] . 

The Wilson line can wind once, twice or more around the periodic temperature direction, 
and the path integral average of the corresponding traces is then written as Trexpifcy, k 
running from 1 to N — 1, and the matric C is supposed to be traceless and diagonal NxN . 

We define now the effective potential U (C) in a box of size L x L y L z as the constrained path 
integral: 



exp-L x L y U(C) = J DAn z , k S(Trexpik^~TrPK(A ))cxp-^S(A) 



(16) 



c We use the phase here as a mathematical device to compute physical obseripgbles like the average of the 
spatial 't Hooft loop. Its physical significance has given rise to a lot of discussion t!J 
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In this definition the area L x L y of the x — y cross section of the box factorizes out. k is 
the winding of the Wilson line. The bar stands for averaging over x and y directions. The z 
dependence in the diagonal matrix C and the Wilsonline is not shown. 

The size of the box is supposed to be much larger than the mass gap of the theory. 

This definition is manifestly gauge invariant, and is invariant under a Z(N) transformation 
of the Aq variable at a given fized timeslice. What does change under such a transformation 
is the value of the Polyakov loop in the constraint. So the argument C is shifted by such a 
transformation, but the value of the effective action stays the same. This invariance is the 
electric Z(N) invariance. 

17(C) has been evaluated up and including two loops Elllil. The calculation is done in 
terms of a saddle point expansion around a background field B in the fourth component of the 
vectorpotential Aq = B + gQa- B is chosen diagonal and the saddle point equations tell that 
B = C. C has N diagonal elements given by C\, C2, • • • , Cm and the constraint J^i Q = 0. It 
lives in the Cartan subspace of the Lie algebra of SU(N). 

To this two loop order the potential is of the form: 



U(C) = _ J2 K{C l3 )(d z C t] f + ]T V{C ij: C ih ) (17) 

i,j ijk 

where CV, = ^ T 3 . Note that these differences are precisely the diagonal elements of the 
adjoint representation of C. So the dependence of the effective action is solely on the adjoint 
representation C a dj of C 

Often we will refer to all of the kinetic term as K(C). 

In this section we will use the effective action to compute the interface tension between 
any pair of Z(N) vacua, with a charge difference A A; = ±1 or what amounts to the same: the 
tension in the area law of the 't Hooft loopa V\ (L) . 

The tension of the 't Hooft loop is computed from eq.(^) by the Gibbs trace: 

TrVx (L) exp - J ~ exp —a\A(L) (18) 



Reexpressing this Gibbs trace in terms of a pathintegral yields precisely the same potential, 

N - 



with boundary conditions C = and C = ^S-TY imposed as discontinuity on a surface spanned 



by the loop L (see fig.(|l|)). 

To see that the average of the loop implies naturally the Wilson line effective action, we 
rewrite the Gibbs trace as a path integral: 

<Vi(L)>= J DAiDA exp{~ Jj dt J d 3 x (d A^ - {D % A ) a - 5(<)<) 2 + (B a ) 2 } (19) 

where we dropped the normalizing factor. The source term af is the coefficient of the linear 
term in E in the exponent of V\(L) in the Gibbs trace. It sits therefore on the surface S(L) 
of the loop, points in the Y direction of the Cartan algebra, and is orthogonal to the surface 
S(L) taken to be in the x-y plane. So it reads on this surface: 

<(x) = S aY 8 iz ^d z e(z) (20) 



d Not anymore true in the presence of dynamical quarks 
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It adds to A the jump as expected for the average of the 't Hooft loop! This jump is 
gauge invariant as we argued below eq.(^J). So we have to find the gauge invariant formed out 
of Aq to meet this invariant boundary condition on Aq. The obvious candidate is P(Aq). 

So the road is clear: introduce the unit into the pathintegral eq.(|l^) by the well-known 
trick: 

1 = J DC(z)J{C(z))n z>k 6(Trexpik^- - TrP k (A Q )) (21) 

with J(C) being a Jacobian. Do the integrations over the potentials in the presence of the 
constraint, as in eq.(16). The pathintegral now reduces to an integral over the profiles C(z), 
with the integrand J(C(z)) exp — L x L y U(C). Letting L x L y to infinity tells us to minimize the 
exponent, with boundary conditions as in fig.(|l]). The phase C jumps at the surface, but its 
gradient is continuous. The gradient is the effective electric field appearing in the kinetic energy 
term of eq.@| 



C/T 





-2it/N 
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Figure 1: Profile of the Wilson line phase for 't Hooft loop in the x-y plane at z = (continuous line). This 
profile minimizes U(C). The broken line is the translate of the right hand branch by the ^-function in eq.(EC) in 
the text, and gives the profile one gets by minimizing U(C) between two Z(N) minima at C = and C = ¥fTY. 



Minimization of the exponent U(C) = K{C) + V(C) means that E{C) = K{C) - V(C) is 
constant in space. Since K = V = 0&tz = ±oo we have E = 0. 

We have drawn the rolling of the system as a function of z in fig.(||). The continuous curve 
corresponds to the continuous curve in fig.(|l]). So this is the profile for the loop. 

In the sections ahead we will compute the surface tension of the 't Hooft loop as follows. 

Once one knows U(C) the same tension follows by choosing the appropriate boundary 
conditions on C (in the case above: C = and C = j^-TY) and minimizing the potential. The 
latter involves finding the path in the Cartan space, that gives an absolute minimum. This 

e The resulting profile looks similar to that of an electrostatic potential due to a dipole layer. There is an 
obvious question here: can one define a phase for the Dirac monopole line, and use it as the magnctostatic 
potential for the profile of the Wilsonloop? 
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-V(C/T) 
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Figure 2: Minimizing the effective action U = K + V with two different boundary conditions. The continuous 
curve corresponds to that of the jump of the 't Hooft loop between 5. and — -j-y . The continuous-dotted curve 
corresponds to the boundary conditions of the domain wall. 
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will be the method followed below and is given in fig. ([!]) by the combination of continuous and 
broken line. The broken line follows by shifting the right part of the profile over This 
undoes the jump at the loop (see eq.(20)). Now the boundary conditions become C — at 
z = — oo and C = ^jfTY at z = oo. One recovers the continuous-dotted curve in fig.(Q) and 
ng.®. 

So the minimum value of U(C) is identical with the two types of boundary conditions. 
This is so because the electric term contains the derivative of the profile, so is the same in both 
cases. The potential has the Z(N) periodicity and hence is the same in both cases. In other 
words: surface tensions of Z{N) domain walls and 't Hooft loops are identical. 



3.1 Effective potential and Wilsonline correlators 

The effective potential can be used as well for the calculation of the correlation of Wilsonlincs 
separated in the z-dircction, and averaged over the x-y directions. This can be seen easily seen 
by Laplace transforming our effective action U{C) with a source term L x L y J.C. It follows 
immediately from the definition of U(C), eq. Jl6|) , that we get the free energy F(J), with the 
source J coupling to the Wilson line. 

exp-L x L y (U(C) - J.C) = J DAexp ( - S(A) + J.P(A Q )) = cxp —L x L y F(J) (22) 

with dcU(C) — J = 0. The dot in J.C means integration over the variable z, and in J.P(Aq) 
integration over all of 3-space. 
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This is nothing but the generating functional of Wilsonline correlators. 
This is important because it permits us to connect the asymptotic behaviour of the profile 
C in the surface tension of the loop to that of the correlator of the imaginary part of two 
Wilsonlines, which is a gauge invariant and non-perturbative definition of the Debye mas£x 
The Wilsonline correlator behaves according to that definition like exp — j-, and so will the 
profile we find by minimization as described above. Ie is the inverse of the Debye mass mu- 
lt will permit us to absorb entirely the 0(g 3 ) term of the surface tension (see section!^). 



4 One loop ratio 

In this section we will compute the surfacetension between two vacua with a Z(N) charge 
difference k. 

From the previous section we know that knowledge of the potential U(C) is not enough to 
compute the tension. We need the right boundary conditions as well. That is done in the first 
part of this section. We will search for the domain wall boundary conditions. 

The potential term lives in the space of traceless real diagonal NxN matrices C, i.e. the 
Cartan subalgebra of SU(N). It has symmetries due to charge conjugation and Z(N) invariance. 
Apart from those there is the periodicity in the C-variables (because the potential is defined 
through the trace of expik^, k — 1, • • • , N — 1) which becomes translation invariance over 2tt 
in the space of the N — 1 variables Q spanning the Cartan subalgebra. 

These invariances lead to the notion of an elementary cell in the Cartan subalgebra. It is 
the smallest simplex in the Cartan subspace that contains all the knowledge about the effective 
potential. _ 

From SU(2) and and SU(3)li-3 one can readily generalize what this cell is. It boils down to 
finding suitable generalizations af-jthe hypercharge Y. There are many equivalent definitions but 
a convenient one is the followingll3: take the diagonal Cartan matrix y\ = Y = (1, 1, • • • 1, —N + 
1) and cyclically permute the diagonal elements to get the Cartan matrices y% = (1, • • • 1, — N + 
1,1), et cetera till un = {—N + 1, 1, • • • , 1). This set obeys the constraint J2k V k ~ 0> an< ^ spans 
all of the Cartan subalgebra {]. 

Obviously 

/ ,2tt . 2tt 
exp(i<— yk) = expil— (23) 

is a centergroup element, whatever the yk- We omitted the unit matrix on the righthand side. 
The corners of the cell are formed by the diagonal matrices 

k 

Y k = ^2 yi = (k,---,k,-N + k,---,-N + k) (24) 
i=i 

with N — k entries k and k entries — N + k. 

Because of eq.(p3|) the corresponding SU(N) group elements expi^Yfe = expik^- are 
centergroup elements. Within this cell there is a hyperplane that reflects all into their 
counterparts Y/v-fc, i-e. represents charge conjugation invariance of the potential term with 
respect to this reflection. Because of Z(N) invariance there are N of these hyperplanes in the 
cell reflecting the arbitrariness of one reference vacuum to another. Only when one breaks 
the Z(N) invariance by introducing fundamental multiplets in the theory this plane of charge 



f It is known that the minima of the potential remain fixed in these points to all orders of perturbation 
theory 
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reflection is fixed as above. So strictly speaking there is a redundancy in our cell due to charge 
conjugation. The reader can find an example of the SU (41 cell in section ^, fig.(^J). 
Now the explicit form of the one loop potential termliJ is: 

V = J dz-TWJ^Biidj) (25) 

y 

B4 is a Bernoulli polynomial(see appendix A). 

To one loop order we have for the kinetic term simply the classical expression for the electric 
field strength: 

K= J dz^Tr{d z C f (26) 

Let us now derive our result (||) for the one loop approximation. We have to minimize 
the sum of K and V going from Yq to Y k . We do this under the assumption that the path of 
minimal action in Cartan space is the rectilinear one: 

C(q) ee Y ok (q) = fa (27) 

where q lies in the interval [0, 1] and parametrizes the path. So on this path CV,( 9 ) = ±q or 0. 

Let us look at the one loop potential V along this path. Now the argument of all Bernoulli 
functions is either ±q or 0. B^q) is even and vanishes at q = 0. We have to sum over all i,j 
combinations that give a non zero result for B^(q) in eq. J25|) . So it is just a matter of counting. 

From eq.(p4[) we see: if i < j we can pick i in N — k ways and j in k ways to get CV, = q. 
Similar for i > j. 

So along the path Y"o/c (q) w e have 

V(C = Y 0k (q)) = k ^~^ V(C = Y 01 (q)) (28) 

for the one loop potential. 

We now evaluate the kinetic term -i^(d z C) 2 , which can be rewritten, using the tracelessness 

of C as 2^1vE 4J (^Cy) 2 - 

Substitute C = Yo k (q), then again there is no contribution from Cy, or it gives the contri- 
bution (d z q) 2 . Exactly the same counting as for the potential terms shows that the sum over i 
and j gives a dependence on k identical to that in the potential term: k(N — k). As it is the 
only dependence on k we have for the ratio for the effective action along the path from to k 
and and along to 1: 

U(C = Y ok (q)) = H ^~^ U(C = Y 01 (q)) (29) 

Note that the ratio is already true for the potential. A fortiori it will be true for the tension. 

The tension o~ k is obtained by minimizing the effective action U along the path Yo k (q). In 
particular for k = 1: 

U = 47T 2 T 2 (JV - 1) J dz(-L-(d z q) 2 + ^q 2 (l - q) 2 ) (30) 



with the well known result I 



u\ = ^ElS (31) 
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The profiles q as function of z are identical for the multiply charged loops. This follows from 
minimizing U{C) and using the result eq.(p9|). This will be the case for two loops as well. Of 
course we would expect the Debye masses to be the same for any number of loops and any of 
the paths Y k(q). 

The counting of contributions leading to the result eq. (|2^) is of course nothing but counting 
the number of non-zero entries in the adjoint representation Yk, a dj(q) as underneath eq. (p7|) . 
Up to a sign they are all equal. 

We emphasize that we took for granted that the straight path realizes the absolute mini- 
mum. It certainly is a local minimum, and only in ££7(4) we are able by numerical means to 
establish the truth of our assumption. 

Now we turn to the two loop ratio. 



5 Two loop ratio 

In this section we show that also the two loop contribution along the straight path eq.(|27|) to 
<Jk has the dependence k(N — k). 

First we analyze the potential term, then the kinetic term. 



5.1 The potential term in two loop order 

We start by recalling the result for the two loop approximation tp_the effective potential V. 
The two loop result is given by the sum of two contributionslHllj : 

Vf = g 2 Y,\f bca \ 2 (B2(Cc)B 2 (C a )-Bl(0)) 



b.c.a 



V P = ig N 



^B 3 (Cij)Bi(Cij) 

i<j 



(32) 



The £?-functions in this expression are the Bernoulli polynomials given in the appendix. 

Vf is the contribution from the free energy graphs, and the constraint in the definition 
eq. (|l^) tells us there is a renormalization of the Wilson line phase Ag. 2 NB\ which adds through 
the lowest order V in eq.(p5|), or rather its derivative B 3 the quantity Vp to the effective action. 
This is useful to keep in mind when we analyze the effective electrostatic action Cs in section ^. 

Let us first look at the contribution of this renormalization of the Wilson line. It has per se 
the same form as the lowest order contribution namely a sum over the index pair ij. So along 
the straight path Yofc this summation gives, as in the lowest order case a factor k(N — k) and 
no other k dependence. 

The contribution of the two loop free energy graphs needs a little more analysis. 

What we have to do is to write out the adjoint indices a, b, c on the f-symbol in terms of the 
index pair symbols ij,jl, li from the Cartan basis. In that basis is a set of N — 1 diagonal 
and orthogonal matrices with norm i and the remaining N(N — 1) Ajy = -^SikSji ■ This is of 
course nothing else than computing with graphs in the double line representation. 

We rewrite in this basis the contribution as: 



V f = . 9 2 ^(5 2 2 (C y ) + 2B 2 (C y )S 2 (0)) 
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+ y( X + J2 )(-B2(Cy) J B 2 (C , i 0+5permutatioiis) (33) 

i<j<l Kj<i 

The first term gives the effect of graphs with one diagonal gluon. It has the same structure 
as the one loop, and hence the now familiar k dependence. 

The second and third term differ by the order of the indices i,j and I in the summation. 
Charge conjugation changes A 1 - 5 into A 3 * and so the first term is the charge conjugate of the 
second. The arguments of the Bernoulli functions only change sign but Bernoulli functions 
of even order are even under sign change of their argument. So along the straight path the 
contribution of the first term is proportional to (iV — k)k{k — 1), because the index i must be 
in the first N — k entries of Yofc(g) = diag(fc, • • • , fc, — N + k, ■ ■ ■ , —N + k), the index j in the k 
last entries of Yofc(<z), so for I there remain only k — 1 entries, because of the ordering. 

The second term is according to the same reasoning proportional to k(N — k)(N — k — 1). 

So adding the first two terms give a proportionality constant k(N — k)(N — 2). Adding all 
three terms of the potential term gives simply: 

V = k(N - fc)^ 2 T 4 B 4 (g)(l - 5^) (34) 

and shows that along the straight path the potential term is just multiplicatively renormalized 
by the two loop terms. This generalizes what was noticed long ago a for the straight path (the 
"q- valley") from the k = to k = 1 minimum, Yoi(<?)- 
So far for the potential term. 

5.2 The kinetic term to two loop order 

The kinetic term one gets by doing a gradient expansion!!. It reads: 

K{G) = 2akt £ + &Y (i,{Cij) + ij(1 ~ Cij)+lE + k^ {dzCij)2 (35) 
ij 

The sum excludes i = j. The ip are the logarithmic derivatives of Eulers T function, and 
develop a pole — for small Cy. They are defined on the interval [0,1] and periodically 
extended. So the sum of the two "0 functions is even under CV, <-> — Cy . Note that for small 
profile Cij is negative. 

For the moment we will ignore this unphysical behaviour for small Cij, and postpone its 
discussion till the section [| We will simply assume in this subsection that 

^(C l j)(d z C % jf = if i,j such that C l3 = (36) 

and justify that in the next section. 

Then, again it is the now familiar counting: along the straight path Y"ofc(<z) w e find easily 
that the kinetic term behaves like: 

K(C = Y Qk ) = k M-^K{C = Y m ) (37) 

Our final task is to minimize the full effective potential U, the sum of the kinetic term 
eq.(k37n and the potential term eq.(|34|) over all possible profiles q(z) along the straight path. 
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This is standard procedure. We take the full potential and write it as a pure square and a 

rest: 

U = J dz(Ki -vi) 2 + 2 J dz(KV)i (38) 

The rest term is independent of whatever profile one takes, since dzK? — 2g \ N K(q)^dq, 
so the rest integral is a definite integral. 
So the minimum of U is reached when 

#3 = yh for all z (39) 

The rest integral then gives the k dependence of the tension as quoted above in eq.(^J). It 
is important that in the rest integral the behaviour of V at q = cancels out the pole in K. 

We have gone into some detail because of an important issue: the minimization goes wrong 
for q — 0(g 2 ), since K becomes then negative. When computed correctly this does not happen 
as explained in the next section. Instead of 0(g 2 ) the correction to the kinetic energy will turn 
out to be 0(g). 



6 The kinetic energy at small values of the profile 

This section gives the justification of ojir procedure in the previous section. This justification is 
based on the effective action approaches. We found it very convenient to use this effective action 



approach for the effective potential of the Wilsonline, the constrained path integral eq.(|16|). 

As a byproduct we compute the well-known first leading correction to the Debye massES. 
It shows the expediency of the method. We also show how the 0(g 3 ) correction to the tension 
in eq.(||) does vanish. The latter is a new result. n 

So the next order that may contribute is 0(g 4 ). Recent lattice datatl confirm that the two 
loop approximation for the surface tension works quite well for reasonably low temperatures. 
Of course only computation of the 0(g 4 ) term will decide, and this is presently being donclia. 

So let us now come to a quantitative analysis of the kinetic term at small values of the 
profile. 

The kinetic term, eq. (|35|) , has been computed in a gradient expansion Since \k\ = mo, 
the Debye mass, the expansion becomes questionable when C ~ gT or smaller. In fact, simple 
powercountingcl tells that for C = 0(g 2 ) the naive expansion breaks down completely already 
at one loop, and this is reflected by the pole ~ g 2 ^ fl. 

6.1 The electric QCD action with a gauge invariant infrared cut-off 



What we .will find by integrating out the hard modes in eq.(16) is an adaptation of a well known 
approachO. 

The adaption is precisely this: in the phase C of the Wilsonline we have a gauge invariant 
infrared cut-off of the propagators. 

So we start from eq.(|l6|), pick a background £ gauge, and integrate out only the hard modes. 
To lowest order one gets the following result: 

C E = Tr(DQo) 2 + ]T TrF 2 nn + -Tr{d.Q + i£{C, Q }) 2 + V(Q , Q; C) + U(C) (40) 

m.n 



9 lt was analyzed some time ago how to remedy the kinetic tern 
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All fields are static and scaled by \/T: = ^ M ^L °- ) . D = d + igz [Q„ F mn is the three 

dimensional gaugefield strength in terms of Q, and 53 = gVT- The gauge condition is the 
static remnant of the original 4d £ background gauge. U{C) is the result we get for U(C) 
incorporating only the hard modes]]. To the order we are interested in the potential term is: 

V(Q , Q; C) = Trm 2 E (C)(Q ) 2 + Tr(i[C, Q]) 2 + 4iTrd z C[Qz,Qo] (41) 

The higher order vertices are of no concern for us in the calculation of the kinetic_cnergy 
for small profile. The potential term has the form of the familiar electric QCD actionO, with 

some important differences however. The propagators of the charged gauge particles Q in 
this 3d theory still have a mass proportional to C%j . Note that this mass is gauge invariant 
and hence assures of a gauge invariant infrared regulator. The Higgs mass matrix m E becomes 
only in the limit of small C equal to the ^g 2 T 2 . There are also additional couplings due to 
the presence of the gauge invariant constraint. These are necessary for the masses to be gauge 
invariant. 

What we need to know is the mass matrix m E (C) for finite values of C. It has three 
contributions. For sake of simplicity we will limit the discussion to the SU{2) case, where the 
matrix C is parametrized as C — diag(7rTg, —nTq). 

The first one is the familiar one that one gets from summing the hard modes with vertices 
and propagators coming obtained from the action with background field C: 

mlM = | g 2 r 2 (l + a3 ,i29 2 ) (42) 

for the diagonal (m 3 ) and off-diagonal (m^) masses. What is important is analyticity in q 2 , 
because we summed only over the hard modes. 

The second and third are due to the presence of the constraint. In what follows we explain 
how they come about. We Fourier analyze the constraint with variables A(z). That adds in the 
exponent in eq. jl^ ) in addition to the action —j^S: 

i J d z^.(±TrP(A )(z) - ilYexpz^) (43) 

We have to perform a steepest descent in both A = A c / + gX qu and the field variables 
Aq = B + gQo and A = gQ. After expansion of the constraint and the action one finds that 
the saddle point is 



TL\ r sin irq 



(44) 

tr 

B = C (45) 



When we develop the constraint one finds: 



-TrP(Ao) - -Trexpz— = --sm{wq) 



JL^(QZ(p0 = 0))2 CO8{7rq ) 



2 



^9lipQltfl M7rq)+om (46) 
2T (p° + 2-KqT) V H> W ' y ' 



h This term is the analogue of the unit operator term of refi 
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The Matsubara frequency is po — ^ttTuq, no integer. For notational reasons we have 
suppressed the z dependence in all the fields. Note there arc no heavy diagonal modes, and 
that the coefficients of the constant modes are different when q ^ 0. 

So there will be a contribution to iriE(q) 2 from the constraint alone by combining the value 
for A c ; in eq. (^) and the expansion of the constraint. In the limit of small q it reads for the 
masses rn%(q) and rn\ 2 (<l) of respectively diagonal and off-diagonal fields Qq the same: 



and for general q: 



rn 



ml(q) = ^ (47) 



2 / x ™;q 



III 



c.3 



2 

c,12 



(q) = — ^cosfa) (48) 
sm(7rg) 

(q) = ^ (49) 



The third contribution comes from the integration over the A^,, variables. 
Their effect is to couple the terms quadratic and higher in Qo from the constraint in eq. 
(Eq) to the conventional vertices present in the action S. 




Qo Qo Q( 



Figure 3: The two diagrams with the mixed constraint/conventional vertices, giving the mass m^ c in eq. (|50|). 
In the loops are only hard modes 

For the mass term TrQ\ these couplings are shown in fig.(0) \ ■ 

The loops contain the hard modes only and give rise to the derivative of the result for the 
one loop potential, eq.(|25|), but with the cubic term omitted (it is due to the zero Matsubara 
frequency). As a result we have for the mass of the diagonal mode: 

<3 = -\ 9 2 T 2 -^- cos(^)(l + 2 9 2 ) (50) 
3 sm(7rq) 

and for the off-diagonal mode: 

m 2 cc , 12 = ~lg 2 T 2 (l + 2q 2 ) (51) 

1 These mixed couplings are important for the consistency of the approach as the following example may 
illustrate. 

There is a contribution from the mixed vertices with the static Qo replaced by the heavy Qo modes with 
p° 0, the last term in eq. When integrated over the heavy modes this gives the diagrams in fig. ^ but 

now with the two lower lines connected to form a loop as well. The upper loops will give again the derivative of 

the one loop potential ~ B^{q) and the loop with the heavy Qo modes gives Bi(q). The latter factor carries the 
gauge dependence necessary to kill theigaLUge dependence in the two loop free energy graphs, like in eq. (p^).This 
mechanism was first noted by BelyaevEj. 
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They are renormalized by the same amount through the hard modes. 
Notice the — sign! And that as expected the integration over the hard modes gives a result 
analytic in C . The final result for the diagonal and off-diagonal masses in the electrostatic 



Lagrangian is obtained by adding eq.([42[), eq. fl49|), eq.(|50[) and eq.(|5l|). The only reason the 
diagonal and off-diagonal masses are different is through the introduction of the constraint. 

The sum of the three types of mass terms is the input in the calculation of the kinetic 
energy in the next section. Summaring this subsection we have for both diagonal and off- 
diagonal masses: 

m%(q) = ^ + 0(q 2 T 2 ) + 0(g 2 T 2 q 2 ) (52) 



in the limit of small q. For the off-diagonal mass 0(q 2 T 2 ) is absent, eq.([49]). 

For small values of q we know from subsection ( ]3 . 1[ ) and from the explicit solution of the 
equation of motion for q, cq.(|39|), that 2t£ = m 2 D , the physical Debye mass. 



6.2 Calculation of the kinetic energy term for small profile with Ce 

This has to be done with the effective Lagrangian eq.(ffOj). The £ gauge choice is a convenient 
choice because the Qo propagators decouple from the Q propagators. We calculate as in the 
4d case the two point function of Gy. Linde's argument u shows that with a magnetic mass 
q = 0(g 2 ) infinitely many diagrams contribute to order g in the kinetic energy with a coefficient 
Ki, in eq.([53|). There is one diagram, shown in fig.(^) that is leading by a log 1 factor. Its 
calculation will be described below. 




Figure 4: Leading contribution to the kinetic energy of the profile (wavy line). Curly line is the Q z , straight 
line the Qo propagator. 



The incoming momentum is k = (0, 0, k z ), the masses in the propagators are rriE(q) for the 
off-diagonal Qo propagator and 2ixTq for the spatial gluon propagator. The vertices are taken 
from the last in term term in V in the electric QCD Lagrangian. It contributes for £ = 1 the 
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term proportional to the logarithm in the kinetic term| : 

(2ttT) 2 ^ , ^ 2 -9 2 T , ik + 2nqT + m E (q) . . 

Here we took specifically g = 0(g 2 ). This is as low as we can go with q in perturbation 
theory. There are many multi- loop diagrams that give all 0(g) contributions embodied in the 
coefficient K\ in eq.(p3[) with this parametric form for q. 

The result of fig. (0) is analytic in k = and has branchcuts on the imaginary axis, starting 
at ±i(m E (q) + 2irqT). Of course we are interested in the straight paths iofc(g) and for SU(2) 
there is only one, k — 1. 

For small q we have q(z) — exp— TOd|z|, as we argued at the end of section^. This means 
that ik — mi). 

The mass m 2 E (q) equals m 2 D plus terms of order 0(g 2 q 2 ) 7 according to eq.(H^). 
So the onset of the branchcuts is at mrj + 2nTq, neglecting 0(g 2 q 2 ). Then the self-energy 
graph, with the external legs at mo, becomes 

g 2 N T { 2mD + 2nTq, 

— log ( — — ) (54) 

2tt mrj 2nl q 

So the final result for the kinetic energy for small values of q is: 

K = ^ ~ {Kl9 + log ^ {dzqf (55) 

Our minimization equation, Ki — V? , tells us that the first term in the kinetic energy is 
balanced by the order one term Vo in the potential term V = Vq + g 2 V2- This gives for the 

Debye mass the well known value tud = y^dT an d that d z — 0(g). The correction to this 
value comes from the second and third term in cq.(|55|). But no term in V can balance these 
terms. So the kinetic term alone must be zero. That is, the 0(g) term has to be absorbed by 
the Debye mass : 

mo = y| 5 T(l + g(±Kt + log ~) + 0{g 2 )) (56) 

as follows from the minimization equation. 



6.3 Why there are no 0(g 3 ) contribution to the surface tension 

An immediate consequence of the last subsection was that the profile realizing the minimization 
equation is such that its kinetic energy to 0(g) is zero. 

So there is no 0(g 3 ) term in the tension 2 j dq(KV)i coming from combining the 0(g 2 ) 
term in VlJLk is absorbed by the correction to the Debye mass! Eq.(|5^) is the well known result 
of RebhanEa for N = 2. 

There is still another possibility for having an 0(g 3 ) contribution: combining O(l) in the 
kinetic energy with 0(g 3 ) in the potential term V. This can be quickly discarded. The reason 
is that a term g 3 Vs(q) is only present in a region where q = 0(g) or less. So integrating over q 
in this region gives an 0(g A ) contribution to the tension. 

J We find that the ultimate logarithmic term is £ independent 
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mo + 2nTq 



Figure 5: The complex k plane and the onset of the branchcuts for the graph in fig. i 



The reader may have noted that in the usual treatment of the log correction to the Debye 

mass the branch cut starts at \J\gT + K, where K is an ad-hoc cutoff of 0(g 2 ). The same 

result as in eq.([56|) is then obtained. But then for small enough g the corrected mass will hit 
the branch cut! In our approach this does not happen. 

Our approach can be easily generalized to the case of SU(N). The diagram in fig.([l]) 
contains as incoming line the field in one of the directions. Some masses of the spatial gauge 
fields will be zero, but it is easy to see, that only the massive ones do couple, with mass 2iiTq. 
This is all we need to find the kinetic energy and the log correction to the Debye mass. 



7 A numerical exploration of the absolute minimum for SU(4) 

That the straight paths are local minima is easy to check analytically. But for the global nature 
of the minima we reverted to numerical search. 

The group ST/ (4) admits 4 minima, Two, Y\ and Y3, are charge conjugate. Their tension 
o 1 has been computed in and is given in our eq.( with N = 4. 

So here we will be concerned with the tension 02 one gets by minimizing U{C), eq.(^8|) and 
eq.(|3|), on any of the paths connecting Y and Y 2 as indicated in fig.(||). 
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Figure 6: Elementary cell for SU(4). The broken line connecting the corners Yq and Y2 is a typical path on 
which we minimized the effective action U(C). Indicated are the directions of the F and G axes defined in the 
text. 



To see the local minimum for the path 102(9) we parametrized the 3d Cartan space of 
£17(4) with the unit vectors \Y 2 and two others, F = -±^(1,-1,0,0) and G = -^=(0,0,1, 1). So 
any point is given by fla + tF + sG. Substitution in the potential shows that indeed a local 
minimum results for r = s = 0. 

To check the global minimum of the effective action minimized on the straight path we 
selected a set of planes all going through the straight path 102(9) and labelled by an angle 9. 

Then we calculated numerically the minimal action on a set of paths having an angle <f> 
with the straight path from lo to Y% . So <j> = is the straight path itself. 

So in fig.(0) each curve corresponds to a given plane. Because of Z(4) and charge conjuga- 
tion symmetry it is sufficient to look at planes in a quadrant < 9 < j. 

It is clear from fig.(^), that the straight path, located at zero slope Taruf) = 0, and present 
in all planes, so in all curves labelled by 9, has globally the minimal action. 

In fig. (Q) we once more plotted the result for the 9 = plane. Note that at the endpoint, 

where the path goes through the point Yi we find 2<7i, and as expected from eq.(0) the value 

3 

2 ■ 

8 A simple and less simplistic approach to the ratios 

We come back to our model of gluonic screened charges we described in section (||) ■ Is it capable 
to reproduce the simple ratio eq.(||)? The first question is: what is the 't Hooft loop with k 
charges? Naively one would say: 

47T f -> 

V k (L) = expi — k / dS.TrEY (57) 
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Tanip 



Figure 7: The minimized action it on a set of paths in a given plane labeled by 8. A path is selected by the 
angle <j> as in fig. (ph. On the horizontal axis Tan<f> is plotted. Endpoints of each curve are where the point D 
hits the Y3Y1 axis in fig. (H). 



just by starting from k singly charged loops far apart and coalescing them. However, this is 
wrong. It is wrong because its corresponding pathintegral expression eq.([l9|) tells us to go k 
steps in the Y direction. Minimizing would then give ka\. But we know already that going 

(7 1 . For a concrete example 



from C 



to ^TYfc gives us a smaller tension a k = k 



(SU(4)) look at fig. g. 

The right answer is therefore: 



(JV-l) 



V k (L) 



47T 

exp* — 



dS.TrEY k 



(58) 



5(L) 



The reader can check that this definition satisfies the fundamental commutation relation eq.(|3|). 
When worked out in terms of the usual pathintegral expression, as in section ||, it tells us to 
follow the minimal path from C — to C — ^-TY k in the elementary cell. 

To see how the ratio comes about we assume that the gluons are screened, free and of 
course in the adjoint representation of the colour group. So we will denote the off diagonal 
gluons by an index pair ij. The diagonal gluons play no role. 



To understand the counting in the case of Vk(L) we use as in subsection 2.1 the adjoint 
Wilsonloop TrPexpi J L ,dl.A a dj looping Vfe(L). As the reader can check with the canonical 
commutation relations: 



V k {L)Wadj(L')V k (Ly 



2tt 

exp {i—Y ktadj )W a dj{L') 



(59) 
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Figure 8: The ratio for paths in the plane = 0. The end point of the curve corresponds to a = 2a\ as 
explained in the text. 

where the hypercharge matrix in the adjoint representation is given by its diagonal elements. 
These are computed from the differences of the fundamental hypercharge matrix Yj. and are 
therefore or N. So the adjoint Wilsonloop commutes with the 't Hooft loop as it should. But 
the pointlike adjoint charges send only half of the flux through the 't Hooft loop, so produce a 
minus sign. And the structure of Yh,adj tells us precisely how to count. 

Again, a gluon given by an index pair ij will have charge or ±A, with respect to the 
U(l) group defined by Y&. Let us call this charge lij(k)(±N). For example, if k — 1, then all 
hj(k) = 0, except when one of the indices equals N, in which case (1) = 1. A gluon with 
generic index pair ij contributes (— to the average of 14 in eq.(p8[). 

The slab of thickness 21 e around the minimal surface spanned by the loop contains gluons, 
with a Poisson distribution P(hj) for each of the off-diagonal ones, Uj runs from to oo. Their 
mean number present in the slab be kj = I, assuming they are all the same. This should be the 
case because of global SU(N) symmetry. Even if the distributions were not Poissonian their 
widths should all be the same for the same reason. Now we find for the average of the loop 
Vfc(L): 

<V k (L)> = UijPikjX-) 1 ^ (60) 
= exp-4fc(A- k)l (61) 



There is a factor 4 because each index pair ij has a conjugate ji which contributes the same 
to the exponent, namely 21 or 0. All what goes in eq.(61) is that a gluon ij with hj(k) = 
contributes a factor 1 to the average, and a factor exp (—2/) if hj(k) = 1. The rest is counting 
of the index pairs that have hj(k) — 1, and follows the reasoning just above eq. (Bp) . 
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A non-Poissonian width will contribute a common (ij independent) factor in front of the 
exponent. 

The result for the ratio eq. @ follows immediately, since I = ti(T)AIe where n(T) = riij (T) 
is the density of a charged gluon ij. 

So also the ratio is a consequence of the gluons being screened, free and in the adjoint 
multiplet of SU(N). 

9 Spatial Wilsonloops 

About the Wilsonloops there is no analytic knowledge, unlike the 't Hooft loops. What we 
know is from dimensional arguments, and from simulationscJ. 

The dimensional argument says that very high temperatures the physics of our system 
is determined by 3d gauge theory with coupling g| = g 2 T . So the surface tension p will be 
proportional to this mass scale squared: 

p~ 5 4 T 2 (62) 

Another important fact from simulations is that in the hadronic phase the surface tension 
is constant from T = to T — T c . This is consistent with the idea that in this phase the cause 
of the area law is colour magnetic flux in the T = groundstate, and that the hadronic gas at 
T ^ does not contribute any such flux. 

Above the critical temperature, according to the dimensional argument and corroborated by 
simulations E-3, the string tension behaves as if the plasma phase itself provides colour magnetic 
flux. 

So the question becomes: what is providing this flux? 

There are basically two answers: the first is macroscopically long flux loops. Their presence 
in the plasma is unlikely since it is difficult to excite them thermally. Only will those loops 
contribute that maybe already present in the groundstate of cold QCD. So they will not cause 
the rise in the Wilson surface tension with temperature. 

The second option is monopolcs. Their fluxes are screened in the plasma. The Dirac 
monopoles with flux are too singular. 

Then there are 't Hooft-Polyakov type monopoles with screened flux. In what follows we 
will concentrate on them. What we have to say here is speculative. 

We will assume there are screened free monopoles at very high T in_the plasma. They 
are supposed to have charge _ An old argument by Goddard et al. suggests that the 
magnetic group is given by SU(N)/Z(N) if the gauge group is SU(N), or SU(N) if the gauge 
group is SU(N)/Z(N). We suppose the monopoles to be in the adjoint representation. They 
are consistent with the presence of free quarks with charge -ft. For N — 2 they have the same 
charge as the 't Hooft-Polyakov monopoles. 

With this hypothesis we can deduce the area law for the Wilson loop, and the ratios for 
Wilsonloops with centergroup charge k. Not surprisingly it will turn out to be identical to 
that of the 't Hooft loops, eq.(^). Not surprising because our input is an adjoint multiplet of 
monopoles dual to the multiplet of gluons. And the counting of the effective gluon charges, as 
in eq. (js^) , has a dual as we will see shortly. 

We start with the area law in SU(2). 

We need a representation of the Wilson loop that clearly shows how it captures the flux of 
an 't Hooft-Polyakov monopole. .— . 

In a recent paper by one of the authorsc3 the following representation of the spatial Wilson- 
loop in pure SU(2) was found by precisely a procedure of continuation of the 't Hooft-Polyakov 
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monopolc flux through a loop L from the Higgs phase to the symmetric phase (the VEV of the 
Higgs goes to zero) and then letting the Higgs mass go to infinity, which is essential to give the 
expression a well-defined meaning. What one finds is that the fundamental Wilsonloop equals: 

W X {L) = /n CTiTiQ Z37i Q ((7,T)exp^ / dS.(B a n a + e abc n a Dn b f\ Dn c ) (63) 

J 2 JS{L) 

= J Dflexpi J dsTry (nA„(s)n _1 + -fld s n) (64) 
The unit vector ft is defined as: 

ft = TrfTVsQ^ (65) 

The integration over ft shows the fluctuations over the directions of the Higgs field that 
come in naturally when letting the VEV of the Higgs go to zero. These representations were 
first derived by Diakonov and PetrovEH. They prove extremely useful in the present context. 
In the first representation we have a clear connection with the integrated flux of an 't Hooft- 
Polyakov monopole. That of the second equation is clearly a rewriting of the Wilson loop. The 
path integral in the first expression is over the unit vector ft defined on any choice of surface 
with as border L. The second path integral is over all SU(2) gauge transformations periodic 
on the loop L. 

In order to compute the average of the Wilsonloop let us assume that our monopoles, being 
screened, are free. This makes only sense if their mean distance exceeds the screening length 

im {]• 

With one monopole present within a distance Im of the loop the first representation gives 
a factor —1 for the loop. 

How does this happen? Like for the electric charges in the presence of the 't Hooft loop 
the monopole is carrying a magnetic charge vis a vis the t%. So the monopole is off-diagonal 
and picks up the difference of the eigenvalues of T3. Note the factor i in front of the exponent. 
The outcome is therefore n for the flux in the exponent. 

Then with the now familiar Poisson distribution P(l) for the monopoles in the slab of 
thickness 21m we find the area law: 

< Wi{L) >=^2(-) l P(l) = exp-2l (66) 

l 

and the mean number of monopoles in the slab is I = 21mAum(T). 

As mentioned in the introduction the density of monopoles is therefore related to the surface 



tension of the Wilson loop and using eq. (62) 



n M (T)~g b T d (67) 
So the density of the monopoles compared to that of the gluons is down by a factor (- — tttt) 3 . 

B\ /J,rp ) 

There is an important consistency check on our assumption that the monopoles are free 
and screened. As we said before, their density should be small enough with respect to the 
screening length. So we should compute the number of monopoles within the screening radius, 
I m iim(T). From the relation between the tension and the density it follows: 

l 3 M n M (T) ~ l 2 M p(T) (68) 



fc See the discussion section for discussion of this point 
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Only for Sll(2) both the screening mass j from the correlation of heavy Dirac monopoles and 
the tensionc3 are numerically available: 

l-J = 0Mg 2 (T)T (69) 
p(T) = (0.369 ± 0.014) V^ 2 (70) 

So we find for SU(2) that the number of monopoles inside the screening radius is ~ 0.19, 
or the mean distance between two monopoles is about twice the screening length 0. 
Finally the pressure will get a positive contribution to the coefficient of 0(g 6 ). 

9.1 Ratios of Wilson surface tensions 

Let us turn to the ratios of surface tensions for SU(N), with N > 3. For that we need a formula 
for the multiply charged Wilsonloop Wk(L). We take again the matrix Y k and write: 

W k (L) = J n a , T , a Dn a k (a, r) exp ± J o!S.(B a n a k + f abc n a k Dn b k A Dn c k ) (71) 

Dn(s)exp— j dsTrY k (n(s)A s (s)n- 1 {s) + -0(s)5 s fi _1 (s)) (72) 
TV J i 

The field n a describes the coset SU(N)/ SU(k)xSU(N - fc)xZJ(l) which has 2k(N - k) 
dimensions. As in the case of the 't Hooft loop it is defined by 

n% = TrY.n^-n- 1 (73) 

where fi is any regular SU(N) gaugetransformation. 

The reader can verify that the fundamental commutation relations |^ are satisfied, using 
the identity TrY k Yi = N 2 mm(k, I) - klN . 

The length of the vector n k (x) does not depend on x: 

a 

The counting of the monopoles that do contribute to the tension p k goes in full analogy 
with that under eq. (p9|) . The difference is that we now ask what the commutation relation is 
between an 't Hooft loop in the adjoint representation V a dj(L) and the Wilsonloop W k (L'): 

V adj (L)W k (L')V-^ (L) = exV i—Y k , adj W k (L') (75) 

with Y k a( ij the Y k charge in the adjoint representation. 

So the loops do commute. But now the pointlike monopoles send only a flux through 
the Wilson loop, so will contribute — 1. 

It is now familiar from section ^| how to proceed with the counting of monopoles that do 
indeed contribute —1. An off-diagonal monopole ij will contribute only when ij is one of the 
2k(N — k) combinations admitted by Y k . The Wilson surface tension becomes proportional to: 

p k - k{N - k)T (76) 



^Obtained in raf. f^ l for a lattice of temporal extent 4 at /3 = 2.82. 
""We neglected in the data a discrepancy in subtraction points. At large enough T this is irrelevant. 
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where I is the average number of ij type monopoles in the slab of width 21m, taken to be the 
same for all. 

We then find the scaling for the ratio: 

HN - k) 

Pk = N-l Pl (77) 



10 Discussion 



In this paper we have looked at the physics of the plasma through the 't Hooft loop and the 
Wilson loop. The former is on a firm basis, the latter of a more speculative nature, but both 
give testable results. For the former we showed: 

i) that the surface tension associated to the 't Hooft loop has a calculable perturbative expansion 

ii) that it obeys, including two loop order and up to order 0(g ) a simple dependence on the 
Z{N) charge in Cartan space of the loop minimized along the straight path 

iii) that the Debye mass absorbs all g 3 effects of the surface tension. 

iv) that the Feynman rules coming from the constrained pathintegral have a natural infrared 
cut-off built in. 



For the lowest order the result ii) was evident, for the two loop result it is due to the 
factorization into two Bernoulli functions. For three loops the factorization is gone, so one has 
to check by different methods wether the simple dependence still holds. We should point out 
that there are statements in the litterature E3 that already the second order contribution has 
non perturbative contributions. That this is not the case is demonstrated in section ||. 

Let us remember that this same simple dependence came out of the plasma viewed as a 
gas of free screened gluons. From that point of view the interactions could easily change the 
behaviour. 

For SU(4) our result for the 't Hooft loop is on a firm basis. We checked that the straight 
path is the global minimum. What is painfully lacking is some analytic understanding of this 
fact in SU(4), and for higher N. What singles out the straight paths for all values of N is that 
the Wilson line effective potential renormalizes multiplicatively along these paths, something 
that was known up to now only for the path from to Y\ giving the singly charged 't Hooft 
loop tension. The ratios at the critical point (which we could obviously not compute with 
the techniques at hand) will give valuable information on the 3D spin model describing the 
universality class of the transition. 

Point iv) suggest doing multi-loop calculations for the free energy with these Feyman rules, 
in particular the perturbative logarithm at sixth order. 

Let us now discuss the more speculative aspects for the Wilson loops. 

For the Wilson loop we could understand on the basis of lattice results that at high T a 
second component of screened and free monopoles can explain the tension and the numbers 
available for SU(2) suggested this is a consistent picture. 

The magnetic contribution to the sixth order pressure has. now a definite sign and is related 
to the surface tension of the Wilsonloop measured by Karschu through Im, a quantity measured 
from the correlation of the Dirac monopole loopslia. 
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The reader might be shocked by the statement that within the screening radius there are 
few monopoles. For the statistical Debye screening for the gluons one has of necessity many 
gluons within the screening radius Ie- I e Ue(T) ~ g~ 3 - 

On the other hand the magnetic .screening has not only statistical origins, because also 
at low temperatures it is still present 113. So we do not see any flagrant inconsistency in the 
assumption that there is a dilute gas of screened monopoles present in the plasma at very high 
T. It is dilute because it only has to explain the "small" Wilson loop tension p(T) ~ g 4 T 2 . On 
the other hand the density of gluons has to explain the "large" semiclassical result for the 't 
Hooft loop tension, a{T) ~ ±T 3 . 

If it is indeed the dominant contribution to the sixth order term in the pressure, it would 
because of its positive sign push the pressure up from the known fifth order, and nearer to the 
Stefan-Boltzman limit. We do not pretend that it sheds light on the perturbatiyeiConvergence of 
the pressure in the plasma phase of SU(N) gauge theory! In fact a recent letter E3 to understand 
the convergence needs a negative sixth order coefficient. This result is not in contradiction with 
ours! The sixth order coefficient consists of the non perturbative contribution and a perturbative 
one. The latter consists of a term coming from the T scale, and a log term with argument Id/Im- 
Knowledge of the perturbative part would clearly be helpful. 

For the ratios of Wilson loops we speculated that these monopoles were in the adjoint 
representation of a magnetic group (SU(N) or SU(N)/Z(N)). Then the same argument that 
worked for gluons gave the dependence on the centergroup charge k of the Wilsonloop. 

In the ratios the unknown constant from the statistical distribution drops out and they are 
the same as those for the 't Hooft loops. 

In order to get a more quantitative grip one should have a precise lattice measurement 
of Im, the sixth order coefficient of the pressure, the perturbative logarithmic term therein, 
and the surface tension of the spatial Wilsonloop. The SU(2) measurement of the magnetic 
screening mass should be done for SU(3), so that together with the spatial Wilson loop results 
one can check the consistency of the monopole model in that case as well. 

A measurement in SU(4) and SU(b) of ratios at high T would settle some of the issues. In 
the case of the spatial Wilson tension at very high T we have the option of simulating in the 
T = 3D theoryfi 

A remark applying to our results for both types of loops: their surface tensions do not 
depend on the representation, only on the N-allity. This follows from our discussion of their 
effect on physical states in subsection 2.1. The surface effects can only come from the 't Hooft 
commutation relations. The latter are determined by the N-allity of the loops. As expected, 
Teper and LuciniLj find the sameihr temporal Wilson loops in T = SU{A) and 577(5) theory. 
This issue is still debated thoughts. 

What happens when dynamical quarks are present in our plasma? As we noted, they are 
consistent with the monopole component through the Dirac condition. So the Wilson loops will 
not be affected by their presence. But the 't Hooft loops will be affected. They will depend on 
the surface S(L), not anymore only on the loop L. These effects can and should be studied. So 
our plasma would have three screened components at high temperatures: gluons, quarks and 
monopoles. The latter are by far the largest screened objects, and in density down by a factor 
(log-^-)~ 3 . At lower T they start to coalesce. 

The results described in ref.0 have bearing on the magnetic activity in the cold phase and 
its behaviour in the hot phase, i.e. on the temperature independent part of the Wilson loop. 
Our monopoles have to do with the rise of the tension of the Wilsonlop. The fact that our 
monopoles are screened may be due to this magnetic activity. 

As a last question: what happens to the ratio in 3d gauge theories? There the answer 
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is: the same as in 4d gauge theories for the analoque of the 't Hooft loop. That is: take the 
correlation between two vortices with strength k and compare it to that with strength k = 1. 
Both have stringy behaviour and the ratio of the stringtensions is as in eq. (Q) . 

For the Wilson loops the situation is different. The tension is caused by the 't Hooft Z(N) 
vortices. If one assumes that vortices with different charge have the same density at high 
enough T we find the ratio is one, but we expect in general a dependence. 

When finishing this paper we were informed about the-.existence of string/brane inspired 
computations of ratios of Wilson loops at zero temperatureES. There the ratio is in the notation 
adopted here: 

pu sin k , , 

V = (78) 

pi sin 

This has in common the liiiearity in present in our result. This result is consistent with 
recent findings on the lattice t£. Of course the physics in both cases is very different so it is 
not clear wether comparison makes much sense. 
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Appendix A. Bernoulli polynomials 

The hatted Bernoulli polynomials appearing in the text are related to the Bernoulli polynomials 
in Ryzhik-Gradshtein, as given below. 

B i ( x ) = ^T*(B i (x) + ^) A.1 

B 3 (x) = yT 3 B 3 (x) A.1 
T 2 



2 

Bx(x) = -— Bx{x) AA 

47T 



with 

B 4 (x) = x 2 (l - \x\) 2 

3 1 

B 3 (x) =x 3 - -x 2 e(x) + -x 

B 2 {x)=x 2 -\x\ 
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=x - -e(x) 

All £>'s are defined mod 1. If the subscript is even, the B is even; if odd the B is odd. So 
in particular B\ is a sawtooth function. 
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